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ABSTRACT 

We present calculations of the dynamics of highly ionized gas clouds that are confined 
by external pressure, and are ionized by AGN continuum. We focus on the gas that 
is seen in absorption in the X-ray spectrum of many AGN and show that such gas 
can reach hydrostatic equilibrium under various conditions. The principal conclusion 
is that the clouds can be accelerated to high velocities by the central X-ray source. 
The dynamical problem can be reduced to the calculation of a single parameter, the 
average force multiplier, (M). The typical value of (M) is ~ 10 suggesting that radi- 
ation pressure acceleration by X-rays is efficient for £/-^Eddington 0.1. The terminal 
velocity scales with the escape velocity at the base of the flow and can exceed it by a 
large factor. The typical velocity for a HIG flow that originates at R — 10 17 cm in a 
source with L x = 10 44 erg s^ 1 is ~ 1000 km s^ 1 , i.e. similar to the velocities observed 
in several X-ray and UV absorption systems. 

Highly ionized AGN clouds are driven mainly by bound-free absorption and 
bound-bound processes are less important unless the lines are significantly broad- 
ened or the column density is very small. Pressure laws that result in constant or 
outward decreasing ionization parameters are most effective in accelerating the flow. 
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1 INTRODUCTION 

Highly ionized gas (HIG) is common in both types of ac- 
tive galactic nuclei (AGN). It is seen as strong absorption 
features in many type-I (Seyfert 1) galaxies (Reynolds 1997; 
George et al. 1998, Kaspi et al. 2000 and references therein) 
and as strong, large equivalent width emission lines in many 
type-II (Seyfert 2) AGN (e.g. Turner et al. 1997; Netzer, 
Turner & George 1998). The situation regarding the high 
luminosity AGN (the bright quasars, see George et al. 2000 
for review and references) is still unclear, mostly because of 
the limited signal-to-noise ASCA observations used in such 
studies. 

The strongest observed HIG features in type-I AGN are 
several absorption edges, at around 1 keV, mostly due to 
Ovn (0.74keV) & Ovm (0.87keV). The deduced hydrogen 
column density is in the range of 10 21 - 10 24 cm , the 
covering fraction is close to unity and the ionization param- 
eter about two orders of magnitude larger than the broad 
line region (BLR) ionization parameter (George et al. 1998) . 
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Recent Chandra observations reveal the presence of numer- 
ous absorption and emission lines (e.g., Kaspi et al. 2000). 
The absorption lines are narrow (~ 150 km s~ )and exhibit 
outflow velocities of ^ 1000 km s" 1 . Currently, there is no 
clear indication regarding the location and hence the mass 
of this component. The very few observations of variable ab- 
sorption features (e.g., Guainazzi et al. 1996; and references 
therein) suggest a dimension not too different from that of 
the BLR. If correct, this would mean that the HIG observed 
in type-II AGN is of different origin since its distance from 
the centre must be larger than ~ 10 pc. 

Many type-I AGN show also narrow UV absorption 
features due to CivA1548 and other resonance lines. In 
most cases the lines are blueshifted with respect to the sys- 
temic velocity, implying outflow with typical velocities of 
~ 1000 km s _1 . In several well studied cases the lines are 
resolved, showing intrinsic width of 100 — 300 km s _1 (Cren- 
shaw et al. 1999; Crenshaw & Kraemer 1999; Srianand 2000; 
and references therein) . There have been several suggestions 
linking the origin and dynamics of the UV and X-ray ab- 
sorption features (Mathur et al. 1994; Mathur et al. 1995; 
Murray & Chiang 1995). This issue is still under discussion 
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(e.g. the new observations of 3C288.1 by Hamann, Netzer & 
Shields 2000). 

Radiation pressure driven flows are thought to be com- 
mon in astrophysical environments. Stellar winds can be ac- 
celerated to their observed velocities 1000 km s _1 ) by 
radiation pressure force (e.g., Castor, Abbott & Klein 1975). 
Broad absorption line (BAL) flows in AGN (BALQSO) can 
be driven to sub-relativistic velocities (~ 30, 000 km s _1 ) by 
the same physical mechanism (e.g., Arav, Li & Begelman 
1994; Murray et al. 1995; Proga, Stone & Kallman 2000) 
although other mechanisms have also been proposed (Begel- 
man, de Kool & Sikora 1991). In both cases, the absorption 
line widths are of the order of the outflow velocity which lead 
to a continuous flow geometry. In this case, the physics can 
be described by the Sobolev approximation which is valid 
provided the velocity changes on length scales over which 
the other gas properties (e.g., ionization state, temperature, 
density) do not change considerably. Obviously not all types 
of flow can be driven by radiation pressure forces, among 
these are jets and supernovae envelopes. 

A second, almost orthogonal approach to the study of 
flows is to consider moving "clouds" or condensations (here- 
after, the "cloud model" ) . Such clouds are confined by means 
of external pressure and are likely to be in hydrostatic equi- 
librium. These models have been studied thoroughly in the 
context of broad line region (BLR) dynamics (e.g., Blumen- 
thal & Mathews 1979; Mathews 1986). The clouds are as- 
sumed to be accelerated ballistically, with no appreciable in- 
ner velocity gradients and their spectra consists of narrow, 
blueshifted absorption lines (e.g., Mathews 1975). Cloud 
confinement and stability have been studied, extensively, 
in many works (e.g., Krolik 1979; Mathews 1986; Math- 
ews 1992). Possible confinement mechanisms include ther- 
mal confinement (e.g., Krolik, McKee & Tarter 1981) with 
its well known problems (e.g., Mathews & Ferland 1987), 
and nonthermal confinement, e.g., by magnetic fields (Rees 
1987; Bottorff, Korista & Shlosman 2000). Despite the many 
unresolved issues, the cloud model has had great success in 
accounting for the observed features of the BLR clouds (e.g. 
Netzer 1990), NLR clouds (e.g., Kraemer, Ruiz & Crenshaw 
1998), and BAL flows (e.g. Arav et al. 1994). 

While acceleration by UV radiation, in AGN, has been 
studied extensively, there are only very few works discussing 
the acceleration by X-ray radiation (e.g., Mathews & Capri- 
otti 1985; Mathews & Veilleux 1989, none in great detail). In 
particular, there is no systematic study of the acceleration of 
gas with properties similar to those commonly observed in 
the X-ray spectrum of many AGN (for a heuristic treatment 
of the problem see Reynolds & Fabian 1995, for a magneto- 
hydrodynamical model with comparison to observations see 
Bottorff, Korista &: Shlosman 2000). As shown below, con- 
siderations of the relevant parameter space within the frame- 
work of the cloud model, result in very different dynamics 
for the HIG clouds compared with previous works. 

This paper investigates the acceleration of the HIG 
clouds. The purpose is to study the micro-physics of the 
gas and the radiation pressure force, and to derive the gas 
velocity under various conditions prevailing in type I AGN. 
This will enable to explain observed spectral features in new, 
high resolution X-ray spectra. 

The paper is organized as follows: in §2 we outline the 
main ingredients of the model. Section 3 describes the cal- 



culation procedure, and present the results concerning the 
cloud structure, the radiation pressure force and the gas mo- 
tion, with a following discussion on the model limitations. 
The main conclusions are summarized in §4. 



2 THE MODEL 

2.1 Basic assumptions 

Our purpose is to solve for the structure and the dynam- 
ics of HIG clouds in AGN. Such gas has been detected by 
ASCA and earlier missions but the evidence used here is 
found mainly in the recent Chandra observations such as 
the one reported in Kaspi et al. (2000) for NGC 3783 and in 
Kaastra et al (2000) for NGC 5548. In both cases there is a 
rich absorption spectrum indicating high level of ionization, 
large column density (~ 10 22 cm~ 2 ) and outflow motion 
of <S 1000 km s" 1 . The lines are barely resolved (Kaspi et 
al. 2001) and the typical Doppler width is (~ 300 km s _1 ). 
The observed properties are, therefore, different from those 
found in stellar winds and in BAL quasars. They can be in- 
terpreted as either due to cloud motion or continuous wind. 
This paper addresses the first possibility and the wind solu- 
tion will be discussed in a future publication. 

The main assumptions regarding the cloud model are 
(see also §3.5): 

(i) The general AGN environment can support various 
type of clouds. Known examples are the BLR and the NLR 
clouds (e.g. Netzer 1990; see a recent study by Kaiser et al. 
2000). Such entities with higher level of ionization cannot 
be ruled out by any known AGN property and have, in fact, 
been addressed, by several recent publications (e.g., Contini 
& Viegas 1999; Bottorff et al. 2000). The HIG region is mod- 
eled as a spherical cloud system of large covering factor. The 
justification is that such regions are very common in low lu- 
minosity AGN (e.g. George et al. 1998). We concentrate on 
thin spherical shells although this is not a requirement of 
the model (see §3.5.) 

(ii) We assume that the clouds are confined by ex- 
ternal pressure P xt- This is required since clouds which 
are not self-gravitating evaporate on sound crossing time 
scales (which is typically much shorter than dynamical time 
scales). While the origin of this pressure is not the topic 
of the paper, we note that confinement by a hot gas (e.g. 
Krolik, McKee & Tarter 1981) is problematic (Mathews & 
Ferland 1987) and magnetic pressure is the most likely expla- 
nation (see Rees (1987) for a general formulation, de Kool & 
Begelman (1995) for a discussion of radiatively driven, mag- 
netically confined flows and Sivron & Tsuruta (1993) for a 
discussion of the magnetically confined cloud model). The 
pressure is assumed to be time-independent and its radial 
dependence is given by 

Pcxt oc iT, (1) 

Typically, -10/3 < a < -3/2 (e.g., Netzer 1990). The 
global spherical symmetry is not critical to our discussion 
and is relaxed in §3.5. 

(iii) We restrict the discussion to geometrically thin HIG 
clouds, AR/R <C 1 where R is the location of the illuminated 
surface (see for example Reynolds & Fabian 1995). This is 
based on variability studies which set rough lower limits on 
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the density (~ 10 6 cm -3 , see Otani et al. 1996; Guainazzi 
et al. 1996) and thus, through the known ionization param- 
eter (see below) on the location of the clouds. The external 
pressure difference AP ex t = P cx t(R) — Pcxt(R + AR) satisfies 
AP cx t/P e xt = —aAR/R <C 1. This assumption is checked 
for self-consistency in §3.5. 

(iv) The clouds are assumed to be in photo-ionization 
equilibrium. This is justified as long as the recombination 
time is the shortest time scale. This is true in our param- 
eter range provided the ionizing continuum is not rapidly 
varying (Krolik & Kriss 1995). 

(v) The dynamical problem is treated by assuming radia- 
tion pressure force, gas pressure force, and gravity. The local 
gas acceleration, a(r), at position r inside a cloud whose il- 
luminated face is at a distance R from the centre, is given 
by 



a(r) = a rad (r) - g 



1 dP si 



p dr 



(2) 



where a ra d is radiation pressure acceleration pressure, g the 
gravitational acceleration, p the gas density, and P gas the gas 
pressure. Note that the dependence of g and a ra( j on R has 
been omitted since R is constant across the geometrically 
thin cloud. Furthermore, we assume quasi-hydrostatic cloud 
structure in the cloud's rest frame. As shown by Blumenthal 
& Mathews (1979), and others, this is an adequate descrip- 
tion of BLR and similar clouds. The clouds considered here 
are illuminated by the same radiation field and are roughly 
at the same location, hence the same general arguments ap- 
ply. The clouds are therefore regarded as moving with a 
certain, distance dependent, bulk velocity and acceleration 
((a), i.e., a(r) = (a)), with no significant velocity gradients. 
Given this, we do not have to treat the full hydrodynamic 
problem but assume that the clouds remain stable and do 
not disperse as they are accelerated. The quasi-hydrostatic 
assumption is checked for self-consistency in §3.5. 



2.2 Ionization structure and radiation pressure 
force 

The ionization structure of the cloud is determined by the 
local radiation field. This is usually described in terms of 
the illuminated face (r = 0) ionization parameter, U , 



V 



I"o ( L »l hv ) dv 



(3) 



where L„ is the monochromatic continuum luminosity 
(erg s _1 Hz -1 ), nn the hydrogen number density at r — 
and c the speed of light. Our aim is to study the acceleration 
of highly ionized clouds. We therefore use the X-ray ioniza- 
tion parameter, U x (Netzer 1996), defined by huo —0.1 keV 
and hv\ = 10 keV. The ionization and thermal structure are 
calculated by ION99, the 1999 version of the photo-ionization 
code ION (Netzer 1996, and references therein). The code is 
very detailed regarding X-ray related processes and solves 
for the temperature profile and the radiation pressure force 
distribution inside the cloud. All our calculations apply to 
fully exposed clouds (i.e., no obscuration or shielding). 

The absorption of continuum radiation by bound- 
bound, bound- free, Compton scattering and free- free pro- 
cesses results in a net momentum transfer to the clouds. 
The relative importance of the various processes depends 
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Figure 1. Contribution of various radiative processes to the nor- 
malized radiation pressure force (relative to the Compton radia- 
tion pressure force) for the standard case with line width given by 
"sound- Note the sharp decline of the bound-bound contribution, 
after ~ 1% of the total column, due to optical depth in the lines. 



on the atomic cross sections and the spectral energy dis- 
tribution (SED). The largest cross sections are for bound- 
bound transitions. These occur at energies lower than the 
ionization energy of the relevant ion where, for the assumed 
AGN-type SED (see below) the value of L„ is larger than 
its value at the ionization threshold. However, large cross 
sections imply large optical depths and self shielding. Thus, 
bound-bound absorption contributes the most at the illu- 
minated surface and decays quickly into the cloud (Fig. Q). 
Bound-free processes are the next largest contributers. The 
Compton cross-section is relatively small, and the optical 
depth is negligible in all cases considered here (column den- 
sity < 1.5 x 10 24 cm -2 ). Thus, Compton scattering con- 
tributes equally (per particle) across the cloud. In almost 
all cases (see however §3.1), free-free absorption contributes 
very little to the radiation pressure force (~ 0.01 Comp- 
ton). We neglect internal line radiation pressure (see Elitzur 
& Ferland 1986) which do not contribute significantly to the 
total pressure because of the low excitation temperature of 
the strongest X-ray transitions. 

We use the force multiplier formalism (e.g., Mihalas & 
Mihalas 1999). The force multiplier is defined as the ratio 
between the radiation pressure force due to all processes 
relative to the force due to Compton scattering, 



M(r) 



n e (r)a T L to t 



(4) 



where Xv{ r ) is the frequency dependent total absorption and 
scattering coefficient, n e is the number density of free elec- 
trons, Ltot is the total, unattenuated (r = 0) luminosity, and 
T v (r) is the position dependent optical depth. This gives 



a ra d(0 = 



n e (r)a T Ltot 
4-KR 2 p(r)c 



M(r). 



(5) 



2.3 Model parameters 

Several parameters influence the solution. The SED detr- 
mines the ionization structure and the radiation pressure 
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force over the various energy bands. We have experimented 
with several SEDs typical of type-I AGN and found only 
little differences regarding the structure and dynamics of 
the clouds. We therefore take as our standard and only 
case an SED characterized by a "UV bump",a OI = 1.4, 
0(0.5-1 kev) = 1-06, and a (1 _ 10 keV) = 0.9. 

The internal line widths determine the line optical 
depth and thus the bound-bound contributions to (M). We 
have experimented with three possibilities: a) Thermal line 
width, vt- b) Gaussian-shaped profiles with line widths for 
all ions equal to the thermal hydrogen line width (denoted by 
fsound). c ) Gaussian-shaped, line profiles, characterized by 
their FWHM (typically several hundred km s^ 1 ) and moti- 
vated by the Crenshaw et al. (1999) and Kaspi et al. (2001) 
observations. Mechanisms responsible for such broadening 
have been discussed, extensively, in the literature (see Bot- 
torff & Ferland (2000) for discussion of MHD turbulences, 
including many older references, as well as many papers 
discussing microturbulences and expansion in stellar atmo- 
spheres) . 

The range of ionization parameter and column density 
is determined by the observations (Reynolds & Fabian 1995; 
Reynolds 1997; George et al. 1998; Kaspi et al. 2001). We 
chose 0.02 < U x < 4 and 10 20 < N H < 10 23 cm -2 . Many 
of the results pertain to a "standard model" defined by 
U x = 10" ' 8 , n H = 10 8 cm -3 , N H = 10 22 cm -2 , with line 
width given by vt- We consider only the following solar- like 
composition, H : He : C : N : O : Ne : Mg : Si : S : Ar : Fe = 



2.4 Method of solution 

We follow the general method described in Weymann (1976), 
Mathews & Blumenthal (1977) and Blumenthal & Math- 
ews (1979) where the cloud structure obeys the hydrostatic 
equation 



p(r) K atl (r) - g- (a)} . 



(6) 



dP gas (r) 
dr 

where (a) is the average (bulk) acceleration. This is satis- 
fied for every distance R. The radiation pressure force is 
changing across the cloud due to optical depth effects. This 
results in a non-uniform gas pressure profile which must 
obey the boundary condition P sas (r = 0) = -P ox t(-R) ~ 
Pext(R + AR) = P g as(r = AT?). The exact distribution de- 
pends on the ionization structure, the level population of 
the various ions, and the temperature at every point inside 
the cloud. 

The cloud accelerates ballistically and we make no use 
of the continuity equation. Nevertheless, the cloud column 
density can change due to the external conditions since it 
can accrete mass, lose mass, or expand along the way. In 
this work we treat the ID problem and consider the fol- 
lowing Nh(R)'- 1) constant mass gas shells which subtend a 
constant solid angle, f2 < 47T, i.e. Q.R 2 Nh = const, 2) con- 
stant column density clouds. The latter model is included 
as an example for extreme cases where the column density 
decreases with distance much slower than R~ 2 . For exam- 
ple, spherical constant mass clouds imply N H oc R 2/3a . This 
allows us to study the sensitivity of our results to different 
column densities. In §3.5 we consider several other implica- 
tions of these assumptions. 



To calculate the acceleration (a), we integrate equation 
([]) with the given boundary conditions (i.e., pressure) at 
both cloud surfaces. Given our parameter space and set of 
assumptions, the contribution of the small external pressure 
gradients to the acceleration via AP ex t / (pAP) is negligible 
(~ 10~ 4 <?) and the bulk acceleration is given by 



(a) 



J Afl p(rKad (r)dr 
f* R p( r ) dr 



(7) 



which depends on R via the boundary conditions, the value 
of Nh and via a r ad(>~)- 

Keeping with the definition of the average acceleration, 
we define the average force multiplier, 



(M) 



f* R n e (r)M(r)dr 
J^ R n e {r)dr 



(8) 



where we have made use of the fact that U x is large and 
hydrogen and helium are fully ionized, thus p(r)/n e (r) ~ 
1.18mff- Thus we obtain the hydrostatic equation in the 
force multiplier representation, 



dPg as (r) n e (r)a T L 

tot 



[M(r) - (A/)] . 



(9) 



dr AtvR 2 c 

The cloud dynamics is therefore reduced to the ballis- 
tic case where the cloud is accelerated coherently with the 
following equation of motion 

dv , , m , 1 crrLtot lllr rr,\\ GA/bh 
V dR = {a(R) ' = 



(M(R)) 



(10) 



v " 1.18mn 4ttcP 2 R? 

where v is the velocity of the cloud and ATbh is the mass 
of the central black hole. Thus, (M) is the only parameter 
which depends on the inner structure of the cloud. 

Defining a new variable, /, which is proportional to 
1/A/bh, by 

1 arLtot/c _ Ltot 



1 1.18m H 4ttGMbh 1.18L odd ' 
simplifies the equation of motion to 

dv _ 1 or Ltot 
V dR 



(<M(R)) - i) 



(11) 



(12) 



1.18ms 4ttcP 2 

In the above formalism, the dynamical problem is decou- 
pled from the cloud structure problem (i.e. the structure is 
independent of the cloud velocity). 

The actual calculations start by assuming a uniform 
density profile and by solving for the temperature and ra- 
diation pressure force profiles. We calculate (AT), solve the 
hydrostatic equation and obtain the first estimate of P gas (r). 
Given the first temperature profile, we obtain a new density 
profile which is used, as an input, for the next iteration. This 
procedure is repeated several times until convergence. Tests 
show that the final results depend very little on the initial 
density profile. The iteration stops when the density profile 
does not change by more than 1% between one iteration to 
the next. 

Given the pressure and ionization profile, the solution 
of the equation of motion is straight forward. The initial 
conditions include the distance from the centre at the base 
of the flow, Ro, -Pext(Po), Nh(Ro), the column density as a 
function of R (i.e., constant column clouds or constant mass 
shells), and the AGN parameters (a, Ltot and I). This allows 
us to solve for the cloud structure and for v(R) at each R. 
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3 RESULTS AND DISCUSSION 

This section describes the results of the calculations regard- 
ing the cloud's internal structure, the average force multi- 
plier (M), and the cloud motion. Unless otherwise specified 
the standard model parameters are assumed. 



3.1 Cloud structure 

We allow the ionization parameter to vary in the range 
10 -1 ' 8 — 10 ' 6 assuming all other parameters remain con- 
stant. The contributions of the various absorption and scat- 
tering processes to M axe computed from the ionization 
structure. As seen in Fig. ^, the force multiplier decreases 
with Ux- The reasons are the fewer bound-bound transitions 
at large U x and the lower available driving flux at high en- 
ergies. 

A common property of all models is the increasing den- 
sity into the cloud (Fig. ^ , due to the temperature decrease 
as a function of r and the constant external pressure. The 
most extreme case shown here is for U x = 10~ ' 8 where 
the density increases by a factor ~ 2 above its initial value. 
This is also the ionization parameter where the ionization 
and thermal structure are dominated by H-like and He-like 
oxygen. 

Fig. |^ also shows the internal pressure profile. The rise 
near r = is caused by the pressure gradient that must 
balance the radiation pressure force which is larger at this 
location due to the contribution of bound-bound processes. 
The peak pressure depends on the difference M(r = 0) — (M) 
and the location r where M(r) — (M). For very large values 
of U x , the peak is lower since the gas is highly ionized and 
bound-bound transitions are less important. 

Finally, since the value of the temperature at the illu- 
minated surface, T, is independent of the column density, 
we can approximate it by T ~ 10 5,5 (7°' 5 which is a good 
approximation for the abundances, SED and the range of 
U x considered here. The same functional form, with differ- 
ent coefficients, provides a good fit for other SEDs (see, e.g., 
Komossa & Meerschweinchen 2000). 

Changes in column density affect the pressure and tem- 
perature profile because of the increased opacity. This is 
shown in Fig. ^. For low column densities, and the standard 
Ux, the clouds have almost uniform ionization structure. The 
ionization structure and the radiation pressure force distri- 
bution change considerably for Nh > 10 22 cm -2 since this 
column density corresponds to an optical depth of ~ f in 
the O vil and O VIII bound-free transitions. 

A common feature for all column densities is the fast 
decline of the bound-bound force multiplier, Mbb, beyond 
Nh > 10 17 cm -2 , due to the large line opacity. Deeper layers 
contribute much more to the bound-free force multiplier, 
Mbf. This behavior is column density dependent since, for 
large columns, the increased bound-free opacity causes lower 
energy bands of the continuum to become more significant 
in producing the local radiation pressure (for the typical 
SED considered here, there are many more photons in these 
bands). The increase in Mbf is the reason for the peculiar 



pressure profile in the Nh > 10 cm models. 

We study clouds with a large range of densities (1 < 
nH < 10 14 cm -3 ) at the illuminated surface. As shown in 
Fig. M, the cloud structure is not very sensitive to nn for 
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Figure 2. Density, pressure and force multiplier across the cloud 
for various values of U x . In all cases Nh = 10 22 cm -2 and the 
lines have thermal profiles. The most pronounced difference is 
in M(r), where for lower U x , bound-bound processes contribute 
more near the illuminated surface with a sharper decline towards 
the outer surface. 
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Figure 3. Cloud structure for various Nh- The cloud becomes 
highly non-uniform with increasing Nh and the radiation pres- 
sure force rises towards the leading, non-illuminated edge (see the 
N H = 10 23 cm" 2 curve) 



nnir = 0) < 10 11 cm -3 since the micro-physics remains es- 
sentially unchanged. All changes can therefore be reduced to 
a scale change in AR. Above nn = 10 11 cm -3 , the micro- 
physics changes considerably due to collisional excitation 
and free-free heating. This results in more uniform temper- 
ature and density profiles and hence changes in Mbb- In 
general, the very high density clouds are characterized by 
higher ionization, higher temperature and smaller M. 

The internal velocity spread is dynamically important. 
Wider profiles decrease the line opacity and increase the 
thickness of the zone where M(r) is dominated by bound- 
bound transitions. This shifts the peak of the pressure profile 
deeper into the cloud (Fig. ^|). A large FWHM does little 
to affect the local ionization and temperature balance, as 
expected from purely radiative considerations. However, the 
increased importance of Mbb results, in low column density 
clouds, in a more uniform pressure profile. 
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Figure 4. Density and pressure profiles for different values of 
nu(r = 0). The shape of the density profile is similar in all cases 
up to flu ~ 10 12 cm -3 , beyond which the micro-physics change 



considerably (see the njj 
side). 
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Figure 5. Cloud structure for several FWHM values of oxygen 
lines. The monotonic change in density, pressure, and M(r) is 
the result of the increasing importance of bound-bound transi- 
tions. Note that for FWHM=3000 km s -1 , the bound-bound pro- 
cesses contribute significantly to M(r) even at the leading (non- 
illuminated) surface. 



3.2 The dependence of (At) on U x , Nh, tlh and 
the line width 

The average force multiplier is the primary property deter- 
mining the cloud velocity. It is therefore important to inves- 
tigate its behavior as a function of the cloud properties. 

The dependence of (At) on the ionization parameter, 
U x is shown in Fig. [| Since the clouds are more ionized for 
larger U m , (At) is more sensitive to higher frequency bands, 
where the incident flux is lower. The result is that (At) is a 
decreasing function of U x for all column densities. For low 
column densities (Nh < 10 22 cm -2 ) the change in (At) is 
roughly linear with U x . For larger columns, the structure 
is more complicated (e.g., Fig. |3J) and so is the dependence 
on U x . A crude approximation of the form log (M(U X )) = 
C\ + C 2 \og(U x ), for several column densities and for line 
widths given by u SOU nd, is 




Figure 6. (M) as a function of the ionization parameter, U x , and 
the column density, Nh- A monotonic decline of (M) with U x is 
evident for all Njj . For U x > 10 5 , the cloud is highly ionized 
and Compton scattering dominates the radiation pressure force. 



for U x < 1. The increase in (M) , for N H > 10 22 cm 
result of the changing ionization balance (see text). 
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(13) 



where (At) n is the average force multiplier for Nh = 
10" cm -2 . This approximates (At) within a factor of ~ 2. 

Fig. ^ also shows the dependence of (At) on the column 
density. This is characterized by a decrease followed by an 
increase with a minimum at Nh ~ 10 22 cm -2 . The minimum 
is mostly due to the fact that at this column density, the 
optical depth due to absorption by oxygen is close to unity. 
For Nh < 10 22 cm~ 2 , the clouds are ionized throughout 
(Fig. ^) and the decline in (At) reflects a decrease in the 
available X-ray flux. For Nh > W 22 cm -2 , the ionization 
structure is such that the leading edge is neutral enough to 
absorb a large fraction of the far UV radiation. Since the 
photon flux at those bands is much larger, it results in an 
increasing (At). The changes in (At) are more pronounced 
for low values of U x because of the SED shape. 

The gas density has only a small influence on the cloud 
structure, and on (At), for n H (r = 0) < 10 11 cm -3 (Fig. []). 
For higher densities, the micro-physics changes (see §3.1) 
and the changes in (M) can amount to a factor ~ 2 between 
n H (r = 0) ~ 10 s cm" 3 and n H (r = 0) ~ 10 14 cm" 3 . As 
shown below, for the P ex t oc R~ 2 , Nh = const, case, this 
corresponds to a change of ~ 40% in terminal velocity. Such 
changes can be critical for marginally accelerating clouds. 

The dependence of (At) on the FWHM of oxy- 
gen lines is shown in Fig. ^ where we consider FWHM< 
16,000 km s , similar to BAL flow velocities. We do not 
show cases with u soun d which are roughly equivalent to 
FWHM~ 70 km s _1 . As expected, broader lines produce 
larger force multipliers (~ 6 times larger for FWHM~ 
3000 km s" 1 compared with thermal case). The effect is 
more pronounced for small U x where there are more lines 
available to absorb the larger continuum flux. For U x = 
10" ' 2 there is a decrease in (At) for FWHMJ; 3000 km s" 1 
due to changes in ionization structure which mask the rel- 
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Figure 7. Average force multiplier as a function of njy(r = 0) 
for the standard model with and different U x . For njj < 
10 11 cm -3 , (M) is roughly constant. The changes in (M) at high 
densities are due to the changes in the micro-physics of the gas. 
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Figure 8. < M > as a function of the FWHM of oxygen lines. 
Note the logarithmic scales, and the incline of (M) (for U x < 
10 -0 - 2 ) towards large FWHM values, which characterise BAL 
flows. 



atively small increase in line contributions. For FWHM<; 
16, 000 km s _1 , (M) is relatively constant since the lines are 
optically thin relative to the bound- free processes. 

While the cloud structure and, therefore, the radiation 
pressure force are relatively insensitive to the exact shape 
of the SED (i.e., a x , a ox ) for the same U x , the exact value 
of (M) depends roughly linearly on the X-ray (0.1-lkeV) 
to total luminosity ratio (equation ^J) which in our case is 
~ 0.04. 



3.3 Approximations for (M) 

Fig. ^a illustrates the most important contributions to (M) 
with respect to the cloud dynamics. While the exact value 
depends on the cloud's structure, the bound-free force mul- 
tiplier, Mbf, is, in general, a good approximation to (M). 
It is possible to estimate Mbf by considering oxygen, the 



most abundant metal. Assuming cross sections of the form 
a = ao(E/ Eo)~ s , where Eq is the ionization energy, ao the 
edge cross section, s ~ 2.5 — 3, and an X-ray spectrum with 
an energy power- law index a x . For an ion X with a number 
density n x we find 



M hf = 



_ ( nx \ ( a ° \ 
1 \n e J \ <jt J 



Le(Eo)Eo 
Ltot 



(14) 



where Le(Eo) is the monochromatic flux at the ionization 
edge. For the range of interest, 10" 1 ' 8 < U x < 10" ' 2 , and 
the chosen SED, the main contributions to the force multi- 
plier are Ovu, Ovm, and Hen. Applying the relevant con- 
stants, and assuming Xq V ii + -Xbvm — 1, we obtain the 
following estimate: 

' Ahc II " 



M M ~1.6 + 1.3X vm + 2 



/ AHc II \ 



(15) 



where 10~ 6 ^ Xnen ^> 10 -4 . This Mbf is within a factor ~ 3 
of (M) for line width given by vt- This result as well as the 
general (all elements) Mbf is significantly different from the 
force multiplier calculated by Mathews & Veilleux (1989) 
who considered gas with much lower ionization stages and 
calculated only the contribution due to the carbon edge. For 
the typical U x and SED considered here this gives a value of 
(M) which is a factor of 50(!) smaller than the one calculated 
by us. 



3.4 Cloud dynamics 

Solving the equation of motion (equation |l^) requires the 
value of {M(R)) for every R. The formal solution, assuming 
v(Ro) = 0, is 



v(0 = V c 



1/2 



(16) 



where £ = R/Rq. The characteristic velocity v c is defined as 



1.18m_fi 2ttcRo 



1350 



(ha) 



1/2 



km : 



(17) 



and the initial distance (Ro) in units of 10 17 cm. The dy- 
namical problem is thus reduced to the calculation of a single 
quantity, (M(R)), at each position along the cloud trajec- 
tory. 

We have studied velocity profiles for two generic cases: 
shell-like clouds with constant mass (QR 2 Nh = const) and 
constant column density clouds (Nh(R) = const). The exact 
numerical solutions are shown in Fig. YLu. The diagram illus- 
trates the lower velocities that are obtained for the constant 
column density clouds since these are more massive (note 
that a comparison between the two is not straight forward 
since (M) is a non-monotonic function of Nh, e.g., §3.3). 
For the standard model, the velocities are of the same order 
as the escape velocity at the base of the flow. In particular, 
velocities of 1000-3000 km s _1 are typical of HIG clouds 
originating in or just outside the BLR. Much larger veloci- 
ties are achieved for higher density clouds that originate in- 
side the BLR (related, perhaps, to the inner accretion disc). 



For example, the standard model with nu 



10 1 



and constant column density results in a terminal velocity of 
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~ 20, 000 km s" 1 for line width given by v S outid . Thus, accel- 
eration by X-ray radiation pressure is capable of producing 
BAL-type flows in AGN. This conclusion applies only to the 
ballistic motion of hydrostatically stable clouds. Wind-type 
flows may be different and will be discussed in future works. 

Three cases of external pressure powerlaws were exam- 
ined for constant column density clouds and for constant 
mass shells, all with uh < 10 cm" " 3 (see Fig. |l(J. The 
main results are explained using the constant column den- 
sity cloud case. 

The first case U X (R) = const (a = -2). Here (M(R)) 
is a constant of motion and equation (|l(]) results in 



'•( />'; = \/ 1 - 
where the final (asymptotic) velocity is 



«HnaI = V c 



(M(U X )) - j 



1/2 



(18) 



(19) 



The final velocity of marginal outflows is therefore very sen- 
sitive to the value of I since l/l ~ (M) (Fig. |^b) The accel- 
eration time-scale, defined as the time it takes the cloud to 
reach 90% of its asymptotic velocity, is 



150- 



ffinal 



-R17 years. 



(20) 



In the two other cases, U x decreases with R (a > —2) 
and U x increases with R (a < —2) (see Fig. |l^). In the first 
case, the acceleration due to the radiation pressure force 
decreases as l/£ to some (positive) power for our range of 
a. For £ 1, the clouds become asymptotically neutral and 
the radiation acceleration drops to zero. 

In the second case, U x increases with R, the clouds 
become more ionized as they move out, bound-bound and 
bound-free processes become less efficient and at large £, the 
only driving force is due to Compton scattering. 

The condition for asymptotic outflow (£ ^> 1) is that the 
cloud velocity exceeds the escape velocity at some distance, 



v(R) > v csc (R) 



arL t 



1.18m// 2-n-cRl 



(21) 



Contrary to other wind solutions, in which the final velocity 
is of the order of v eS c(-Ro), the final velocity in our model 
depends on the value of [M(R)l — 1] (equation |l^) which 
can be much larger than unity. 

Fig. [To] also shows that a fair fraction of the final veloc- 
ity is obtained at small £. This is the reason for the similar 
terminal velocity in the a = —2 and the a = —3/2 cases. A 
different behaviour is seen for a — —10/3 where the gravi- 
tational force exceeds the radiation pressure force at £ ~ 3, 
resulting in a marginal outflow. 

The dynamics of constant mass gas shells follows the 
same trend for the various external pressure profiles, but 
with higher terminal velocities. This results from the de- 
crease in column density with distance and the fact that for 
Nh < 10 22 cm -2 , (M) is a monotonically decreasing func- 
tion of N H (see Fig. |). 

So far we have considered radial motion only. Other 
AGN components, such as the BLR and the accretion disc, 
are thought to be rotating around the central mass, and a 
plausible assumption is to associate a rotational motion to 
the HIG. Introducing Keplerian velocities (see also §3.5) at 




(b) 








1=0.2 _ 




1=0.2 




~ - - . J=0.1 




1=0.15 


1 S ~i=u-oa. _ 


j^- 1=0.14 



3500 



-2000 E 
1500 * 
1000 



500 



10 100 1 

!j=R/R n 



10 



100 



Figure 9. Velocity profiles for U x = 10 -0 ' 8 , N H = 
10 22 cm" 2 , R = 10 17 cm, L = 10 45 erg s" 1 , 1 > 1, and P cxt tx 
R~ 2 . Left: v(£) for different approximations of (M): (1) Exact 
solution for thermal line width (i>t)- (2) like (1) but line width 
is given by v Bound , (3) (M) = M(r = 0), (4) (M) = M bf (r = 0). 
As illustrated, the cloud's dynamics is dominated by bound-free 
absorption for both vj< and i) a0 und h ne profiles. Right: v(£) for 
different values of I. line width is given by ■Uaound- Note the large 
differences in i;fi na i for l/l ~ (M). Solid lines are for clouds with 
pure radial motion. Dash-dot lines are for clouds with asimuthal 
(Keplerian) velocities at £ = 1 (see text). Super Eddington lumi- 
nosities are included for completion. 
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Figure 10. The dependence of v(£) on P e xt for a case of 
U x = 10" - 8 , N H (Ro) = 10 22 cm" 2 , R = 10 17 cm, L = 
f0 45 erg s — 1 , 1 = 0.3. Constant Njj cloud profiles in solid lines, 
and constant mass shell (TV/f ex R~ 2 ) in dashed line. The escape 
velocity, v CBC (R), is shown for comparison. For a = —10/3, grav- 
ity takes over due to the rapid expansion of the cloud. Note the 
critical points beyond which the hydrostatic and the thin shell 
assumptions are no longer valid (see text). 



Ro, while conserving angular momentum along its trajec- 
tory, will decrease the effective gravity (in equation [l2], the 
term l/l will be replaced by 1/Z(1 — I/O) thus increasing 
the final velocities (see Fig. ^b). 

The obtained dynamical results differ significantly from 
the few previous calculations involving acceleration of clouds 
by X-ray radiation pressure force. It is an order of magni- 
tude difference from the Mathews & Capriotti (1985) results 
(see their table 1) and this is not surprising since previous 
works did not consider the relevant parameter space for the 
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HIG and, therefore, did not include the relevant radiation 
absorption mechanisms. 

3.5 Limitations and extensions of the model 

We have made two critical assumptions: thin shell clouds 
(which is always true, by construction, at Ro) and hydro- 
static equilibrium at all locations. Below we derive general 
expressions for a critical, normalized distance, £ c , beyond 
which one or both of these assumptions fail. 

The critical distance for the thin-shell approximation is 



oo 



P-a > 1 
otherwise 



(22) 



where we have used nn oc P C xt- Here, (3 — 0, — 2 for constant 
column density clouds and constant mass shells, respectively. 
Since the initial conditions are such that ARo/Ro <C 1, we 
see that in general £< thm ) ^> i. 

A quasi-hydrostatic cloud structure can be obtained 
provided the cloud's sound crossing time is shorter than 
the time it takes the external conditions (i.e., external pres- 
sure) to change considerably (e.g., Blumenthal & Mathews 
1975) as a result of the cloud's motion. This time is roughly 
Ro/vfinai- The quasi-hydrostatic approximation is therefore 
valid provided £ <C £c hyd ' where 



.(hyd) 



/ "sound Rp \ 
V Vfinal ARqJ 



(23) 

"final AR 1 

In cases where /3 — a < 0, the cloud structure and dynam- 
ics can be treated by our model for all £ > 1 provided 
i?o«sound/(A_Ro«finai) > 1. In conclusion, cases for which 



£ c = min(£, 



(thin) t (hyd) 



) < 1 cannot be treated by our model. 



Considering specific examples, we note that for the stan- 
dard model, with L45 = 1, Rn = 1, I = 0.3, U x (Ro) = 



10" 



and q = —2, numerical calculations give £ c 



20. 



At £ = 4 the cloud is already moving at 90% of its final 
velocity. For two other cases of a = -2 (N H oc R- 2 ) and 
a = —1.5, the approximations hold throughout the region. 
The case of a = —10/3 is different due to the rapid expan- 
sion of the cloud. In this case, for Nh = const, our model 
approximations break at small £ (£ c ~ 4). 

Throughout this paper we assumed an external pres- 
sure that depends only on R. This is a natural assumption 
for the constant mass (thin spherical shell) case. For the 
other case, where the column density decrease with distance 
is slower than Nh oc RT 2 , we can imagine a situation where 
the cloud subtends a small solid angle (fi < 4n) and the ex- 
ternal pressure profile is angle dependent. Clouds with finite 
lateral extents are prone to evaporation through their rims 
due to pressure gradients inside the cloud (e.g., Mathews & 
Veilleux 1989, equation 1). This sets a lower limit on the 
cloud's lateral extents and conversly on the spatial angle it 
subtends of 

(^2^) 2 . (24) 





-7— Z W~ 2 
4n 



"final 



Filamentary structure, and different dynamics in different 
directions, can be introduced into our formalism. The treat- 
ment of the general two dimensional problem is beyond the 
scope of this paper (for simple pancake shaped clouds see 
Blumenthal & Mathews 1979). 



In this work we have considered general confinement, 
having in mind external magnetic pressure. Confinement by 
hot thermal gas is another possibility that has been dis- 
cussed, extensively, in the literature. The presence of highly 
ionized, dilute gas between the clouds cannot be ruled out. 
Relative velocities between such gas and the cloud will result 
in drag forces that tend to decrease the final velocities. The 
cloud-gas boundary may be subjected to various instabilities 
which may affect the cloud structure considerably and ulti- 
matly lead to its destruction. This complex situation cannot 
be solved using the methods adopted here. 



4 CONCLUSIONS 

We have made detailed numerical calculations of the dynam- 
ics of highly ionized gas clouds that are in pressure equilib- 
rium with external magnetic pressure and are ionized by a 
typical AGN continuum. The self-consistent calculations in- 
clude the ionization and thermal equilibrium of the cloud, 
the radiative transfer and the one-dimensional hydrostatical 
equilibrium solution. The principal conclusion is that such 
highly ionized clouds can be accelerated to high velocities 
by means of radiation pressure acceleration. The dynamical 
problem can be reduced to the calculation of a single param- 
eter, the average force multiplier, (M). The terminal velocity 
of the flow scales like the escape velocity at its or igin but can 
exceed this velocity by a large factor (~ \J (M)). The cloud 
velocity is sensitive to the confining pressure profile (a), es- 
pecially for marginal flows. In particular, flows that originate 
just outside the BLR in objects with L to t — 10 45 erg s _1 , 
will reach velocities of 1000 — 3000 km s _1 , similar to the 
velocities measured for the X-ray and UV absorption lines. 
Highly ionized AGN clouds are driven mainly by bound-free 
absorption and bound-bound processes are less important 
unless significant non-thermal line broadening or very low 
column densities (< 10 20 cm" 2 ) are involved. 

We thank Nahum Arav for useful discussions, and the 
referee for some valuable comments. We acknowledge finan- 
cial support by the Israel Science Foundation and the Jack 
Adler Chair of Extragalactic Astronomy. 
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